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SECTION  1 


INTRODUCTION 


1.1.  BACKGROUND 

Forced  vibration  tests  were  conducted  on  a model  buried  arch  and  on 
two  buried  rectangular  structures  at  the  site  of  the  ESSEX  Phase  3 Event  near 
Fort  Polk,  Louisiana.  The  tests  on  the  arch  were  the  subject  of  a previous 
theoretical  study  (Reference  1)  which  investigated  how  resonance  was  suppressed 
when  the  structure  was  covered  with  backfill.  A series  of  finite  element  analyses 
showed  that  radiation  damping  was  significant  for  the  prevailing  backfill  condi- 
tions, and  that  the  amount  of  radiation  damping  varied  from  about  10  to  20  per- 
cent of  critical  over  a reasonable  range  of  backfill  wave  speeds.  The  present 
study  considers  the  response  of  the  rectangular  structures  from  a similar  view- 
point. 

The  geometry  of  the  rectangular  structures  and  backfill  and  instrumenta- 
tion layout  are  shown  in  Figures  1-1  and  1-2.  These  structures  were  subjected 
to  forced  vibration  testing  after  the  ESSEX  Phase  3 Event.  The  original  back- 
fill was  excavated,  the  structures  were  inspected  and  found  to  be  undamaged. 

Forced  vibration  tests  were  conducted  first  on  the  uncovered  structures.  Then 
sand  backfill  was  placed  by  raining  from  a height  of  5 feet  and  the  tests  were 
repeated.  The  tests  and  results  are  described  in  detail  in  Reference  2.  Typi- 
cal results  of  the  tests  are  shown  in  Figures  1-3  and  1-4  in  the  form  of  imped- 
ance (F  /v  ) and  quadrature  (F  times  peak  acceleration  times  sine  of  the 
max  max  max 

phase  angle  between  them)  versus  frequency.  When  the  driving  frequency  approaches 
a natural  frequency  of  the  structure,  there  is  a valley  in  the  impedance  plot 
(large  vffl  ) and  a peak  in  the  quadrative  plot  (large  peak  acceleration).  One 
observation  is  that  the  stiffer  structure,  3D  (Jl/h  = 4),  exhibits  resonance 
for  both  the  covered  and  uncovered  conuitions.  In  contrast,  the  covered  arch 
exhibited  no  resonance  at  all  and  the  covered  rectangular  structure,  3B  (£/h  = 10), 
exhibited  weak  resonance.  The  second  surprising  result  is  that  the  resonant  fre- 
quencies of  the  covered  structure  are  almost  identical  to  those  of  the  same  struc- 
ture without  backfill.  This  conflicts  sharply  with  the  results  of  modal  analy- 
sis of  protective  structures  in  soil  which  shows  that,  if  the  soil  is  assumed  to 
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be  bonded  to  the  structure,  the  addition  of  soil  mass  to  the  structure  decreases 
its  embedded  frequencies  relative  to  those  in  vacuum.  Analytic  results  obtained 
from  the  arch  show  this,  but  they  cannot  be  confirmed  because  resonant  frequen- 
cies were  not  observed  for  the  physical  structure. 

Analysis  of  free  vibration  of  the  covered  arch  showed  that  the  amount 
of  radiation  damping  varied  from  about  10  to  20  percent  of  critical,  depending 
on  the  wave  speed  assumed  for  the  backfill.  Since  the  seismic  wave  speed  for 
the  sand  is  reported  as  oeing  30  to  45  percent  of  that  for  the  compacted  native 
backfill,  it  is  anticipated  that  the  radiation  damping  for  the  rectangular  struc- 
tures is  less  than  that  for  the  arch.  This  would  help  to  explain  why  resonance 
was  suppressed  in  the  arch  and  not  in  the  rectangular  structures,  and  hence 
is  one  of  the  main  subjects  of  the  present  study. 

The  observations  that  resonance  appears  in  the  covered  case  of  rectangu- 
lar structure  3D,  that  it  is  muted  in  the  case  of  rectangular  structure  3B  and 
that  it  is  completely  suppressed  in  the  case  of  the  covered  arch  signify  that 
frequency  and  damping  parameters  for  equivalent  single-degree-of-f reedom  (SDOF) 
models  cannot  easily  be  obtained  from  forced  vibration  tests  without  further 
analysis.  The  results  also  point  out  that  deformational  modes  are  strongly  in- 
fluenced by  local  backfill  conditions  and  that  knowledge  of  the  effects  of  back- 
fill is  essential  to  accurate  analytic  modeling. 

1.2.  OBJECTIVE 

One  objective  of  the  present  study  is  to  investigate  the  effect  of  vari- 
ations in  backfill  wave  speed  on  radiation  damping  for  the  rectangular  structure, 
3D  (Jl/h  = 4),  and  to  evaluate  equivalent  fractions  of  critical  viscous  damping. 
Structure  3B  is  much  stiffer  than  the  arch  and  represents  radiation  damping  in 
the  high  frequency  range  of  deformational  modes.  The  second  objective  is  to 
use  the  damping  values  in  an  equivalent  SDOF  model  subjected  to  airblast  from 
a hypothetical  explosion  in  order  to  determine  the  effect  of  including  radiation 
damping  on  the  response. 

1.3.  SCOPE 

A study  was  made  of  the  effect  of  backfill  on  the  response  of  a rectangu- 
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lar  structure  subjected  to  transient  loading  using  three-dimensional,  elastic 
finite  element  simulation  techniques.  Natural  frequencies  and  mode  shapes  of 
the  structure  were  obtained  before  and  after  emplacement  of  backfill.  Dynamic 
step-by-step  analysis  of  transient  loading  designed  to  excite  lower  deforma- 
tional  modes  of  the  covered  structure  was  also  performed. 

As  in  the  case  of  the  arch,  the  chief  hypothesis  was  that  assuming 
stiffer  backfill  has  little  effect  on  the  natural  frequencies  of  deformational 
modes  of  the  structure  and  that  it  increases  the  radiation  damping.  In  the 
case  of  the  arch,  a crude  model  of  friction  between  the  backfill  and  the  struc- 
ture was  included  in  order  to  investigate  the  relative  importance  of  friction 
and  radiation  damping.  It  was  found  that  the  latter  dominates  the  response, 
so  that  friction  was  not  included  in  the  present  study.  In  the  previous  study, 
the  effect  of  hysteretic  material  damping  was  also  considered  and  was  shown  to 
be  small  compared  with  radiation  damping  for  the  low  amplitude  excitation  con- 
sidered here  and  produced  in  the  forced  vibration  tests.  Accordingly,  materi- 
al damping  was  neglected  in  the  present  study;  in  other  cases  where  the  stress 
level  in  the  soil  is  higher  material  damping  should  be  considered. 

Nonlinear  SDOF  models  of  the  roof  slab  for  both  rectangular  structures 
3B  and  3D  were  developed  using  existing  resistance  functions  reported  in  Ref- 
erence 3.  Membrane  effects  wei^  considered,  although  the  threshold  displace- 
ment for  these  effects  was  not  reached  in  the  present  analyses.  The  SDOF 
models  were  amended  to  include  a viscous  damping  term  whose  coefficient  was 
varied  within  the  limits  indicated  by  the  finite  element  analyses  to  determine 
the  effects  of  radiation  damping  on  response  of  the  roof.  Since  there  were 
strong  arching  effects  on  the  roof  occurring  during  a HEST  test  on  structure 
3B,  several  loading  functions  were  assumed  in  order  to  bound  the  possible  range 
of  responses.  These  were  based  on  measurements  of  overpressure  in  the  HEST 
cavity  and  on  measurements  of  interface  stresses.  These  studies  were  performed 
as  nearly  as  possible  as  ordinary  vulnerability  analyses,  except  that  a signi- 
ficant amount  of  viscous  damping  was  assumed.  In  this  way  the  effects  of  in- 
cluding or  omitting  radiation  damping  on  targeting  analysis  are  evaluated. 
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ACCELEROMETER 


STRUCTURE  3B  STRUCTURE  3D 


36-1/2"  26-1/2"  25 


38-1/2"  33-1/4"  31 


BOX  STRUCTURE 


Figure  1-2.  Geometry  and  Instrumentation  for  Arch  and  Rectangular  Structures 
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Figure  1-4.  Impedance  and  Quadrature  for  Rectangular  Structure  3D 
Before  and  After  Placement  of  Backfill  Determined  From 
Tests  with  Vibrators  In  Phase. 
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SECTION  2 


FINITE  ELEMENT  SIMULATION  OF  BURIED  RECTANGULAR  STRUCTURE 

2.1.  INTRODUCTION 

The  first  step  in  investigating  the  dynamic  properties  of  the  soil- 
structure  system  is  to  obtain  natural  frequencies  and  mode  shapes  for  the 
lower  deformational  modes  of  the  structure.  Since  the  true  soil-structure 
system  is  semi-infinite  in  extent,  it  does  not  possess  classical  normal 
modes.  Nevertheless,  a reasonable  approximation  to  local  deformational  modes 
can  be  obtained  by  fixing  the  outer  boundaries  of  the  finite  element  grid. 

The  calculated  natural  frequencies  and  mode  shapes  are  compared  with  measure- 
ments in  order  to  qualify  the  analytic  model. 

Next  a transient  response  analysis  using  direct  integration  methods 
was  performed.  An  energy-absorbing  boundary  condition  of  the  Lysmer-Kuhlemeyer 
type  was  used  at  the  subsurface  boundaries  of  the  soil  island  in  order  to  simu- 
late a semi-infinite  soil  domain.  The  wave  speed  in  the  backfill  was  varied  to 
investigate  the  effects  on  radiation  damping. 

2.2.  NORMAL  MODES  OF  STRUCTURE 

The  geometry  of  the  analytic  model  used  in  the  present  study  is  shown 
in  Figure  2-1.  If  the  access  hatch  in  the  roof  is  neglected,  two  planes  of 
symmetry  exist  in  the  structure,  and  this  is  assumed  in  the  finite  element 
model.  This  elminates  anti-symmetrical  modes  and  frees  computer  capacity  to 
represent  accurately  signals  induced  in  the  soil  by  vibration  of  the  structure. 
Since  the  natural  frequencies  of  the  rectangular  structure  are  considerably 
higher  than  those  of  the  arch  previously  considered,  the  wave  lengths  are 
shorter  and  the  required  discretization  is  finer.  The  finite  element  discreti- 
zation of  the  structure  is  shown  in  Figure  2-2.  Hexahedral  continuum  elements 
are  used  for  the  structure  as  well  as  for  the  soil.  As  shown  in  Appendix  I,  a 
preliminary  modal  analysis  of  the  roof  was  performed  to  determine  optimum  dis- 
cretization; the  same  discretization  was  then  adopted  for  the  floor  and  walls. 


The  properties  of  the  soil  outside  the  backfill  which  are  assumed  in 
the  analysis  are  based  on  initial  tangent  moduli,  Reference  1.  This  is  appropri- 


ate  for  simulating  cases  of  forced  structure  vibration  because  the  strain  ampli- 
tudes of  signals  reaching  the  soil  are  small.  In  the  transient  analyses,  the 
properties  of  the  backfill  were  varied  over  a range  deemed  appropriate  to  rea- 
sonable backfill  conditions.  In  the  modal  analysis  of  the  covered  structure, 
the  backfill  properties  are  assumed  on  the  basis  of  initial  tangent  moduli  mea- 
sured in  laboratory  tests  on  Cook's  Bayou  sand.  The  stress-strain  data  for  the 
compacted  native  soil  backfill  and  for  Cook's  Bayou  sand  are  shown  in  Figure  2-3. 
The  properties  are  summarized  in  Table  2-1. 

The  natural  frequencies  and  mode  shapes  were  obtained  with  a Rayleigh- 
Ritz  procedure  in  which  the  initial  shape  functions  were  a set  of  orthogonal 
harmonic  functions  suggested  by  the  theory  of  plates.  Modes  of  the  uncovered 
structure  were  computed  by  placing  the  structure  in  the  excavation  without  any 
backfill.  Modes  of  the  covered  structure  were  extracted  by  adding  soil  elements 
having  the  backfill  properties  assumed  in  Case  1.  Table  2-2  compares  natural 
frequencies  of  the  analytic  model  with  those  measured  in  the  forced  vibration 
tests.  There  is  good  agreement  only  for  the  237  Hz  (214  Hz  calculated)  mode 
of  the  uncovered  structure.  There  is  fair  agreement  for  the  same  mode  of  the 
covered  structure  (180  Hz  calculated).  It  is  significant  that  the  calculated 
frequency  of  this  mode  is  lower  for  the  covered  than  for  the  uncovered  case. 

This  is  due  to  the  accession  of  soil  mass  to  the  structure  which  more  than 
compensates  for  the  added  stiffness.  However,  the  measurements  show  no  change 
in  frequency  when  backfill  is  added.  The  simplest  explanation  is  that  in  the 
physical  experiment  the  sand  backfill  does  not  follow  the  deflections  of  the 
structure  due  to  arching  and  hence  there  is  little  or  no  accession  of  soil  mass. 
In  general,  the  measured  frequencies  are  higher  than  those  calculated.  Several 
anti-symmetrical  modes  were  also  measured,  which  cannot  be  calculated  with  the 
present  model.  The  mode  shapes  are  illustrated  in  Figures  2-4  and  2-5. 

2.3.  TRANSIENT  ANALYSIS  OF  COVERED  STRUCTURE 

The  mechanism  of  damping  which  is  considered  in  the  present  analysis  is 
radiation  associated  with  wave  propagation  in  a semi-infinite  domain.  A Lysmer- 
Kuhlemeyer  energy-absorbing  boundary  condition  is  imposed  on  subsurface  bounda- 
ries of  the  soil  island  in  order  to  approximate  such  a domain.  About  four  velo- 
city-time cycles  in  the  first  mode  of  vibration  are  complete  before  waves  reach 
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the  boundaries,  at  which  point  most  of  the  signal  is  absorbed.  Damping  is  in- 
vestigated by  observing  decay  in  the  amplitude  of  free  vibration  of  various 
points  on  the  structure.  Motion  is  induced  by  applying  a series  of  line  loads 
to  the  roof,  floor  and  wall  of  the  structure  as  shown  in  Figure  2-6.  This 
distribution  corresponds  roughly  to  the  third  deformat ional  mode  of  the  covered 
structure  shown  in  Figure  2-5c,  whose  calculated  frequency  is  about  180  Hz. 

The  time  history  of  the  applied  loads  is  a triangular  pulse  whose  duration  is 
equal  to  0.0028  second  (about  the  half  period  of  the  180  Hz  mode).  The  total 
duration  of  the  calculations  is  0.055  second,  or  about  20  times  the  duration 
of  the  applied  load.  This  load  distribution  and  time  history  produce  a trans- 
ient response  which  is  strongly  dominated  by  the  180  Hz  mode,  with  the  result 
that  it  is  easier  and  more  meaningful  to  interpret  percent  critical  damping 
from  the  log  decrement  of  displacements  in  successive  cycles  than  if  several 
modes  affect  the  response.  The  load  pattern  is  quite  similar  to  the  Ritz  load 
function  used  to  generate  the  180  Hz  mode  in  the  extraction  of  normal  modes. 

The  backfill  properties  considered  in  the  present  parametric  study 
are  summarized  in  Table  2-3.  Case  1 corresponds  to  the  best  estimate  of  the 
properties  in  the  physical  experiment.  Case  2 corresponds  to  backfill  wave 
speeds  twice  those  in  Case  1. 


2.3.1.  Radiation  Damping 


The  results  of  the  transient  vibration  analyses  are  shown  in  Figures 
2-7  and  2-8  in  terms  of  velocity-time  histories  at  the  midpoints  of  the  roof, 
floor  and  sidewalls.  They  show  that  as  stiffer  backfill  is  assumed,  the  peak 
velocity  decreases  slightly.  A frequency  shift  is  difficult  to  detect,  because 
damping  so  dramatically  reduces  the  late-time  velocity  in  Case  2. 

The  fraction  of  critical  damping  is  evaluated  from  the  log  decrement 
of  peak  displacement  in  successive  cycles  according  to  the  expression: 


where 


2Trr| 

/F-n2 


X , X , = peak  displacements  at  cycle  n,  n + 1,  respectively 
n n+l 

h = fraction  of  critical  damping 
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The  displacement-time  histories  are  shown  in  Figures  2-9  and  2-1C.  Since  the 
transient  response  of  the  model  is  not  purely  harmonic,  it  is  to  be  expected 
that  the  value  of  n will  depend  to  a degree  on  where  on  the  structure  the  peak 
displacements  are  observed  and  for  which  cycle.  In  order  to  obtain  a representa- 
tive value  for  the  fraction  of  critical  damping,  it  is  estimated  from  two  differ- 
ent cycles  and  at  three  locations. 

The  estimated  fractions  of  critical  damping  are  given  in  Table  2-3. 

They  show  that,  for  the  backfill  properties  based  on  the  initial  tangent  moduli 
for  Cook's  Bayou  sand,  damping  is  about  15  percent  of  critical.  This  is  slight- 
ly greater  than  the  damping  (12  percent)  obtained  for  the  buried  arch  based  on 
tangent  moduli  for  native  compacted  backfill.  Table  2-3  also  shows  that,  for 
backfill  properties  based  on  two  times  the  recommended  initial  Cook's  Bayou 
moduli,  the  damping  increases  to  about  22  percent  of  critical.  This  is  also 
more  than  the  comparable  value  for  the  buried  arch  (17  percent). 

The  finding  that  the  equivalent  fraction  of  critical  damping  is  calcula- 
ted to  be  higher  for  the  rectangular  structure  than  for  the  arch  is  surprising 
in  view  of  the  fact  that  the  assumed  backfill  is  stiffer  for  the  arch.  The  find- 
ing also  tends  to  conflict  with  the  experimental  observation  that  resonance  was 
observed  during  forced  vibration  of  the  rectangular  structure  (suggesting  low 
damping)  whereas  resonance  was  suppressed  in  the  arch  (suggesting  higher  damp- 
ing). Other  factors  may  have  played  an  important  role,  as  is  discussed 
below. 

The  values  of  critical  damping  obtained  in  this  analysis  are  consistent 
with  those  obtained  for  a more  flexible,  buried  rectangular  structure,  Refer- 
ence 4.  There,  the  dynamic  load  factor  for  the  response  of  the  structure  to  a 
simulated  airblast  was  calculated.  A single-degree-of-f reedom  model  having 
the  same  frequency  characteristics  was  developed  and  viscous  damping  was  in- 
creased until  the  same  dynamic  load  factor  was  obtained,  which  occurred  at  frac- 
tions of  critical  damping  between  20  percent  and  40  percent.  The  values  obtained 
in  the  present  study  are  at  the  low  end  of  this  range,  which  is  consistent  with 
the  lower  wave  speeds  assumed  here  (550  fps  and  777  fps  versus  1,500  fps  in  Ref- 
erence 4). 


A question  arises  whether  radiation  damping  which  is  measured 
over  several  cycles  is  relevant  to  targeting  analysis  of  protective  struc- 
tures in  which  monotonically  increasing  deflections  are  of  principal  interest. 
Radiation  damping  of  deformational  modes  occurs  wher.  deformations  of  the 
structure  induce  stress  and  velocity  gradients  in  the  adjacent  soil  which  are 
propagated  away  in  the  form  of  waves.  A crucial  element  in  radiation  damping 
is  for  the  structural  members  to  do  work  on  the  soil  so  that  energy  of  vibra- 
tion leaks  out  of  the  structure.  For  the  case  of  airblast  loading  on  the  roof 
of  the  structure,  it  is  clear  that  the  principal  responses  of  the  walls  and 
floor  during  the  first  quarter  cycle  arise  from  deflections  and  rotations  which 
are  compatible  with  the  roof  deformation.  This  statement  needs  to  be  qualified 
in  the  case  of  a high  speed  backfill  where  the  downward  traveling  wave  exerts 
appreciable  stress  on  the  sidewall  during  the  period  of  interest.  We  restrict 
the  discussion  of  the  first  quarter  cycle  of  the  lowest  roof  bending  mode. 
Figure  2-5c,  and  to  low  wave  speed  backfill.  The  soil  does  work  on  the  roof 
and  the  walls  do  work  on  the  soil.  If  the  soil  beneath  the  floor  is  capable 
of  following  its  inward  motion,  then  the  floor  too  does  work  on  the  soil.  A 
SDOF  model  of  the  roof  of  the  type  used  for  targeting  analysis  considers  only 
the  roof  slab,  so  that  all  of  the  work  done  by  the  soil  on  the  roof  remains 
as  kinetic  and  internal  energy  of  the  roof.  Leakage  of  energy  through  the 
walls  and  floor  is  neglected.  Seen  from  this  perspective,  the  first  quarter 
cycle  of  loading  is  as  much  influenced  by  radiation  damping  effects  as  any 
other  cycle.  Furthermore,  a simulation  of  the  type  shown  in  Figure  2-6  induces 
much  the  same  pattern  of  deformations  in  the  soil  and  structure  as  are  induced 
by  airblast  loading  on  the  roof.  The  fact  that  it  requires  several  cycles  to 
obtain  a measure  of  the  energy  leakage  is  merely  a nuisance,  but  the  processes 
of  energy  leakage  through  the  walls  and  floor  are  the  same  in  later  cycles  as 
in  the  first  quarter  cycle.  In  the  finite  element  simulation,  the  extent  that 
leakage  through  the  floor  is  exaggerated  by  the  ability  of  the  interface  in 
the  model  to  support  tension,  the  analysis  exaggerates  radiation  damping.  In 
the  opinion  of  the  authors  the  degree  of  exaggeration  is  small  in  the  cases 
which  have  been  considered  so  far. 

2.3.2.  Arching 

Arching  is  the  phenomenon  whereby  the  stress  exerted  on  a buried 
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flexural  member,  such  as  a roof  slab,  varies  across  the  span.  Arching  modifies 
the  effective  load  on  the  member,  and  hence  it  is  important  in  assessing  vul- 
nerability. Arching  can  occur  due  to  dynamic  interaction  in  elastic  soil-struc- 
ture systems  and  due  to  inelastic  behavior  of  soils,  especially  those  with  a 
large  friction  angle.  Although  the  present  study  was  not  aimed  at  investigating 
arching,  the  geometry  of  the  problem  is  conducive  to  it.  Accordingly,  the  inter- 
face stresses  normal  to  the  roof  along  its  transverse  and  longitudinal  center- 
lines  were  examined  to  see  whether  arching  occurs.  Cases  of  soft  and  stiff  back- 
fill are  considered.  The  results  of  this  examination  are  shown  in  Figures  2-11 
through  2-14. 

The  variation  in  stresses  along  the  transverse  centerline  for  the  soft 
backfill  case  is  shown  in  Figure  2-11.  Here,  the  stresses  are  highest  in  the 
center  of  the  span  and  slightly  lower  toward  the  sidewall.  A similar  trend  is 
observed  for  the  stiff  backfill  case.  Figure  2-12,  in  which  damping  reduces  the 
frequencies  of  oscillation.  In  these  cases  the  transverse  variation  in  stress 
is  negligible. 

The  variations  in  stresses  along  the  longitudinal  centerline  for  the 
soft  and  stiff  backfill  case  are  shown  in  Figures  2-13  and  2-14.  These  figures 
show,  as  do  Figures  2-11  and  2-12,  that  the  stresses  are  larger  at  the  center 
of  the  span.  The  reason  for  this  apparently  is  that  the  interface  stresses  vary 
with  the  velocities  and  displacements  of  the  roof;  since  these  are  larger  at 
midspan  than  at  the  supports,  the  stresses  are  also  larger  at  midspan. 
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Soil  Laver  1 


mt 


Soil  Layer  2 
Concrete 

Sand  Backfill 


Figure  2-1.  Elevation  View  of  Transverse  Section  Through  Finite  Llei.ient 
Model  of  Rectangular  Structure  3D. 


22 


Figure  2-2.  Finite  element  model  of  structure  (for  clarity,  soil  is  not 
shown) . 
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Stress-Strain  Properties  of  Compacted  Native  Backfill  (Arch)  and  Cook's  Bayou 
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Displacement-Time  Histories  at  Various  Points  on  Rectangular  Structure  Subjected  to 
Loading  as  Shown  in  Figure  2-6,  Case  1 (Soft  Backfill). 


Displacement-Time  Histories  at  Various  Points  on  Rectangular  Structure  Subjected  to 
Loading  as  Shown  in  Figure  2-6,  Case  2 (Stiff  Backfill). 


Figure  2-11.  Vertical  Stress-Time  Histories  at  Soil-Structure  Interface  Above  and  Adjacent  to 
Transverse  Centerline  of  Rectangular  Structure,  Case  1 (Soft  Backfill). 


Vertical  Stress-Time  Histories  at  Soil-Structure  Interface  Above  and  Adjacent  to 
Transverse  Centerline  of  Rectangular  Structure,  Case  2 (Stiff  Backfill). 
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Figure  2-13.  Vertical  Stress-Time  Histories  at  Soil-Structure  Interface  Above  and  Adjacent  to 
Longitudinal  Centerline  of  Rectangular  Structure,  Case  1 (Soft  Backfill). 


Figure  2-14.  Vertical  Stress-Time  Histories  at  Soil-Structure  Interface  Above  and  Adjacent  to 
Longitudinal  Centerline  of  Rectangular  Structure,  Case  2 (Stiff  Backfill). 


TABLE  2-1 


MATERIAL  PROPERTIES  USED  IN  TRANSIENT 
ANALYSIS  OF  RECTANGULAR  STRUCTURE 


SOIL 


Depth 

c 

s 

c 

P 

Case 

(ft.) 

fps 

fps 

1 — Initial  Tangent  Moduli 

Layer  1 

0-3.5 

659 

1,296 

Layer  2 

3.5-10.8 

983 

1,953 

Backfill 

0-8.2 

311 

538 

2 — 2X  Initial  Tangent  Moduli 

(Backfill  Only) 

Layer  1 

0-3.5 

659 

1,296 

Layer  2 

3.5-10.8 

983 

1,953 

Backfill 

0-8.2 

622 

1,076 

STRUCTURE 

E = 3.8  x 106  psi 
V = 0.18 

p = 0.000225  lb.  -sec.  z/in.  4 

Roof,  floor  and  sidewall  thickness — 11  inches 
Endwall  thickness — 16  inches 
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TABLE  2-2 


CALCULATED  AND  MEASURED  FREQUENCIES  OF 
RECTANGULAR  STRUCTURE  3D  (l/h  = 4) 

UNCOVERED 


Measured 

Calculated 

Frequency 

Frequency 

Mode 

(Hz) 

Mode 

(Hz) 

1 

112 

1 

237 

2 

214 

3 

235 

2 

310* 

4 

333 

5 

353 

3 

415 

4 

465* 

5 

470 

6 

521 

7 

550* 

8 

660* 

9 

790 

TABLE  2-2  (CONTINUED) 


COVERED 


Measured 

Calculated 

Frequency 

Frequency 

Mode 

(Hz) 

Mode 

(Hz) 

1 

88 

2 

172 

1 

237 

3 

180 

4 

271 

5 

293 

2 

310* 

3 

415 

4 

465* 

5 

470 

*Ant isymmetr ical  Mode — Not  Calculated  by  Model  with  Two  Planes  of  Symmetry 
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SECTION  3 

NONLINEAR,  SINCLE-DEGREE-OF- FREEDOM  MODEL  OF  ROOF  SLAB 

3.1.  INTRODUCTION 

The  studies  described  in  Section  2 provide  estimates  of  damping  para- 
meters to  be  incorporated  in  nonlinear  single-degree-of-f reedom  (SDOF)  models 
of  the  type  used  in  targeting  analyses.  In  order  to  demonstrate  the  signifi- 
cance of  including  reasonable  amounts  of  damping,  an  existing  procedure  for 
analyzing  slabs  by  means  of  SDOF  models  was  adapted  and  parameter  studies  were 
performed  for  an  airblast  loading  from  a HEST  test.  This  section  describes 
the  resistance  functions  and  transformation  relationships  for  the  SDOF  model 
and  the  results  of  the  parameter  study.  Comparison  is  made  between  the  res- 
ponse of  this  model  and  measurements  from  a HEST  test  on  structure  3B  (2./h  = 10). 
Response  of  structure  3D  (£/ h = 4)  to  hypothetical  airblast  loading  is  also  ob- 
tained. 

3.2.  EQUIVALENT  MASS,  STIFFNESS  AND  LOADS  FOR  SDOF  MODEL  OF  SLAB 

The  chief  failure  mechanism  for  shallow-buried  rectangular  structures 
subjected  to  airblast  loading  is  collapse  of  the  roof.  Hence  the  vulnerability 
of  the  structure  is  assessed  by  reducing  it  to  a slab  having  the  same  dimensions 
as  the  physical  roof.  This  slab  is  then  converted  to  an  equivalent  SDOF  model. 

The  static  pressure-deflection  relationship  is  developed  from  yield  line  theory 
for  reinforced  concrete  slabs.  The  assumed  yield  line  patterns  are  shown  in 
Figure  3-1.  A typical  resistance  function  is  illustrated  in  Figure  3-2.  It  is 
divided  into  several  regimes,  depending  on  the  state  of  cracking  in  the  concrete 
and  yielding  in  the  steel.  The  expressions  for  the  various  threshold  levels  of 
equivalent  concentrated  load  and  displacement  constitute  the  static  resistance 
function  for  the  slab.  In  the  present  study,  the  resistance  functions  for  elastic, 
elastic-cracked,  elastic-plastic  and  plastic  regimes  were  adopted  from  Reference  3. 
A correction  for  membrane  effects  was  also  included  when  the  deflection  exceeded 
the  thickness  of  the  slab. 

The  distributed  mass,  loading  and  static  resistance  of  the  slab  are  con- 
verted into  equivalent  discrete  quantities  for  use  in  dynamic  analyses  by  means 
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of  transformation  factors.  These  take  into  account  the  elastic  deformed  shape 
and  the  changes  which  occur  as  plastic  hinges  form.  Transformation  factors 
were  adopted  from  Reference  5.  Favorable  agreement  was  obtained  with  models 
based  on  simpler  transformation  factors  (Reference  6);  hence  it  appears  the 
results  are  insensitive  to  minor  changes  in  these  factors. 

The  equation  of  motion  for  the  nonlinear  SDOF  model  is  written  as 

follows : 

mx  = P + F (3-1) 

where 

m,  F(t)  = equivalent  mass,  external  force,  respectively 

P = internal  resisting  force.  For  undamped  systems,  P 
is  uniquely  a function  of  x,  p(x).  For  damped 
systems  p = p(x)  + cx 

x,  x,  x = displacement,  velocity,  acceleration  of  SDOF  model; 

also  correspond  to  similar  quantities  at  midpoint 
of  span  of  physical  structure 

In  the  computer  program  used  in  the  present  study,  the  equation  of  motion  is 
integrated  using  a fourth  order  Runge-Kutta  technique. 

The  idealized  resistance  functions  overlook  the  effect  on  membrane 
forces  of  earth  pressure  exerted  on  the  walls.  Also,  the  transformation  fac- 
tors for  the  load  do  not  account  for  soil  arching.  Further  unresolved  questions 
are  what  is  the  effective  span,  whether  the  support  conditions  for  the  roof  slab 
should  be  assumed  to  be  fixed  or  simply  supported,  and  whether  shear  failure  is 
adequately  accounted  for  in  the  resistance  functions  and  the  transformation  fac- 
tors. It  is  not  the  purpose  of  the  present  study  to  try  to  overcome  any  of 
these  shortcomings,  but  exclusively  to  study  the  effects  of  including  or  omitting 
damping.  Hence,  the  modeling  approach  described  above  is  adopted;  the  results 
of  parametric  studies  using  various  amounts  of  damping  for  the  roof  slabs  of  rec- 
tangular structures  3B  and  3D  are  reported  below. 

3.3.  NONLINEAR  SDOF  ANALYSES  OF  ROOF  SLABS 

The  effect  of  damping  on  SDOF  models  was  investigated  by  varying  damping 
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ratio  as  a parameter  for  given  loading  conditions  and  comparing  the  resulting 
displacement-time  histories.  Two  loading  conditions  are  considered  which  are 
based  on  measurements  in  the  August  1977  Fort  Polk,  Louisiana,  HEST  test  on 
structure  3B.  One  is  the  pressure-time  history  in  the  HEST  cavity,  shown  in 
Figure  3-3,  and  the  other  is  the  normal  stress  exerted  on  the  roof  at  the  cen- 
ter span,  shown  in  Figure  3-4.  The  former  represents  a reasonable  upper  bound 
on  impulse.  The  latter  includes  effects  of  significant  arching  and  as  such 
represents  a reasonable  lower  bound  on  impulse.  These  loading  functions  are 
applied  to  models  of  both  structures  3B  and  3D,  even  though  only  structure  3B 
has  been  tested.  It  is  probable  the  arching  which  was  observed  in  structure 
3B  would  be  different  in  structure  3D  so  it  may  not  be  typical  of  a field  test 
situation  to  use  the  loads  measured  on  structure  3B  for  analysis  of  structure 
3D.  Further,  stress  on  the  roof  of  structure  3B  probably  varies  in  a way  which 
was  not  detected  because  there  were  only  three  soil-roof  interface  gages.  These 
difficulties  could  be  overcome  with  some  elaborate  models  and  more  measurements, 
but  the  objectives  of  the  present  project  require  that  the  simplest  models  be 
used  and  this  is  what  is  done. 

The  displacement-time  histories  for  the  SDOF  model  of  roof  of  structure 
3B  (clamped  edges)  subjected  to  loading  based  or.  interface  pressure  are  shown 
in  Figure  3-5.  Four  damping  ratios  corresponding  to  0,  10  percent,  15  percent 
and  20  percent  of  critical  were  considered,  which  nearly  covers  the  range  of 
damping  ratios  calculated  in  the  previous  study  on  the  arc!  , Reference  1,  and 
the  present  study  on  structure  3D.  The  damping  coefficient  in  the  SDOF  model 
is  calculated  on  the  basis  of  the  following  expression: 

c = n (2 v'Tcrn)  (3-2) 

where 

k = stiffness  of  the  elastic-plastic  regime  of  the  SDOF  model 

m = mass  of  the  SDOF  model 

0 = fraction  of  critical  damping 

As  Figure  3-5  shows,  there  is  a major  effect  when  the  viscous  damping  equal 
to  10  percent  of  critical  is  included;  the  effect  increases  moderately  when 
the  damping  is  increased  to  20  percent  of  critical.  The  total  effect  for  20  per- 
cent damping  is  to  reduce  maximum  displacement  to  about  40  percent  of  that  in 
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the  undamped  case.  Threshold  values  of  displacement  which  separate  the  various 
regimes  of  behavior  such  as  elastic-cracked,  elastic-plastic  and  fully  plastic 
are  indicated. 

The  computed  displacement-time  history  cannot  be  compared  with  data  be- 
cause no  measurements  of  transient  roof  sag  were  made.  However,  a permanent  de- 
flection or  roof  sag  of  0.5  in.  was  measured  by  surveying  after  the  test  and  it 
is  tempting  to  try  to  compare  this  aspect  of  the  data  with  the  calculation.  There 
is  considerable  uncertainty  in  defining  the  unloading  branch  of  the  resistance 
function.  As  is  shown  in  Figure  3-2,  the  elastic-plastic  and  elastic-cracked 
slopes  are  reasonable  lower  and  upper  bounds,  respectively,  on  the  unloading 
slope;  as  the  figure  suggests,  the  uncertainty  in  these  slopes  can  lead  to  a 
large  uncertainty  in  permanent  roof  deflection.  These  unloading  slopes  were  ap- 
plied in  the  20  percent  damping  case  considered  in  Figure  3-5  (structure  3B,  in- 
terface pressure-time  history  loading).  The  corresponding  resistance-deflection 
path  is  shown  in  Figure  3-6.  If  unloading  is  assumed  to  be  governed  by  the  elas- 
tic-plastic slope,  the  calculated  permanent  displacement  would  be  0.24  in.;  if 
unloading  is  assumed  to  be  governed  by  the  elastic-cracked  slope,  the  calculated 
permanent  displacement  would  be  0.74  in.  These  values  bound  the  measured  value 
of  0.5  in. 

Results  of  the  analysis  of  structure  3B,  with  loading  based  on  pressure 
in  the  HEST  cavity,  are  shown  in  Figure  3-7  for  two  assumed  values  of  critical 
viscous  damping.  There  is  much  more  impulse  in  this  loading  function  than  in  the 
one  based  on  interface  pressure  and  hence  the  calculated  sag  is  much  greater. 

Since  the  absolute  velocities  are  higher  in  the  case  with  HEST  cavity  pressure  in- 
put, the  damping  contribution  to  F,  Equation  (3-1),  is  also  significantly  greater. 
Thus  the  reduction  in  roof  sag  due  to  increasing  damping  from  10  percent  to  20 
percent  of  critical  is  greater  in  this  case  than  in  the  one  based  on  interface 
pressure  loading. 

Similar  results  were  obtained  for  the  model  of  structure  3D  when  it  was 
subjected  to  the  same  loads.  Results  for  the  loading  based  on  interface  stress 
measurements  made  during  a test  of  structure  3B  (and  hence  not  necessarily  appli- 
cable to  structure  3D)  are  shown  in  Figure  3-8.  The  velocities,  displacements 
and  effect  of  including  damping  are  significantly  less  for  structure  3D  than  for 
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structure  3B.  This  is  to  be  expected  since  structure  3D  is  much  stiffer  and 
stronger  than  structure  3B.  As  before,  greater  roof  sag  and  greater  effect  of 
including  damping  is  observed  when  the  REST  cavity  pressure  is  used  as  input, 
Figure  3-9. 

The  interpretation  of  these  findings  is  complicated  by  the  significant 
uncertainties  involved  in  transforming  a physical  structure  into  a SDOF  model. 
These  include  the  effective  span,  the  edge  support  conditions,  the  applicability 
of  yield  line  theory  and  how  to  incorporate  arching  effects.  There  is  also  un- 
certainty concerning  the  fundamental  step  of  representing  radiation  damping  in 
the  physical  structure  by  a constant  percentage  of  critical  viscous  damping 
based  on  the  elastic-plastic  stiffness. 

In  spite  of  these  uncertainties,  the  examples  described  above  indicate 
clearly  that  even  such  light  damping  as  10  percent  of  critical  significantly 
reduces  the  response  calculated  by  SDOF  models. 
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Figure  3-5.  Displacement-Time  History  for  SDOF  Model  (4  Different  Damping 
Ratios)  and  for  Midpoint  of  Roof  for  Structure  3B  (£/h  = 10) 
Subjected  to  Measured  Interface  Stress-Time  History. 
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Figure  3-''.  Displacement-Time  History  for  Structure  3B  (Thin  Box  l/h  = 10) 
Applied  Load  Based  on  HEST  Overpressure. 


DISPLACEMENT  (IN) 


SECTION  4 


SCALING  RULES  FOR  RADIATION  DAMPING  ASSOCIATED 
WITH  VIBRATIONS  OF  A CYLINDRICAL  CAVITY 

The  equations  governing  frequency  and  radiation  damping  characteris- 
tics for  a cylindrical  structure  embedded  in  an  infinite  medium  have  been  formu- 
lated so  that  basic  relationships  among  wave  speed  and  frequency  can  be  under- 
stood for  each  mode.  Such  relationships  provide  insight  into  how  to  scale  or 
extrapolate  knowledge  of  radiation  damping  gained  from  specific  finite  element 
anlyses  to  new  cases  involving  different  geometry  and  material  properties.  Sca- 
ling rules  indicate  how  frequency,  wave  speed  and  size  of  structures  can  be  varied 
to  increase  or  decrease  radiation  damping.  For  example,  if  two  structures  having 
different  frequencies  of  deformational  modes  are  embedded  in  identical  backfill 
and  soil,  it  is  helpful  to  know  how  natural  frequencies  of  the  structures  are 
modified  and  which  is  more  highly  damped.  Moreover,  if  radiation  damping  associ- 
ated with  one  of  the  structures  is  known,  scaling  rules  may  suggest  how  to  esti- 
mate radiation  damping  for  the  other.  Cylindrical  geometry  is  adopted  because 
it  is  the  closest  simple  approximation  to  a plane  section  perpendicular  to  the 
long  axis  of  the  rectangular  and  arch  structures  considered  in  the  present  and 
previous  analyses.  The  work  described  below  assumes  that  the  relationships  which 
apply  to  the  cylinder  also  apply  qualitatively  to  more  complicated  geometry  and 
material  properties  such  as  a rectangular  structure  or  an  arch  embedded  in  back- 
fill with  nonlinear  properties.  Further  field  and  laboratory  testing  and  possi- 
bly computations  are  needed  to  confirm  this  assumption. 

The  following  analysis  involves  coupling  a cylindrical  shell,  represented 
by  its  flexural  modes,  to  a cylindrical  cavity  in  an  infinite  elastic  domain. 

The  coupled  equations  of  motion  are  expressed  in  terms  of  a frequency  equation 
for  each  mode  shape  (rigid  body,  n = 1,  ovaling,  n = 2,  etc.).  The  solution  is 
restricted  to  lower  modes  (n  < 5)  where  coupling  between  flexural  and  extensional 
modes  can  be  neglected,  Reference  7.  The  frequency  equation  has  several  roots 
for  each  mode  shape,  corresponding  to  the  multiple  frequencies  and  damping  for 
modes  of  the  cavity.  The  dominant  mode  can  be  selected  by  evaluating  and  compar- 
ing the  participation  factors  for  each  solution.  The  present  study  carries  the 
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solution  to  the  point  of  examining  the  response  in  each  mode,  including  fre- 
quency characteristics  and  damping.  Synthesis  of  the  modes  could  be  performed 
to  obtain  a complete  dynamic  solution.  Reference  7,  but  this  is  not  our  goal. 
The  overall  solution  to  this  problem,  which  is  given  below,  is  not  new.  A form 
of  the  solution,  using  doubly  asymptotic  methods,  is  available  in  Reference  8 
and  was  used  in  evaluating  the  effects  of  embedment  on  frequencies  and  damping 
ratios  of  a cylinder  in  Reference  1.  However,  Reference  8 is  not  intended  for 
this  purpose;  a more  convenient  solution  is  given  below  in  the  form  of  a chart 
which  is  intended  to  be  used  by  targeting  analysts. 

4.1.  THEORETICAL  DEVELOPMENT 

Consider  a thin-walled  elastic  cylindrical  shell  with  circular  cross 
section,  embedded  in  an  elastic  continuum  as  shown  in  Figure  4-1.  Assume  that 
the  shell  and  continuum  are  in  welded  contact  and  that  a state  of  plane  strain 
exists.  Further  assume  that  any  free  surface  or  material  discontinuity  is  far 
from  the  cylinder  and  that  the  shell  deformation  is  due  primarily  to  in-plane 
flexure  (in  contrast  to  membrane  extension).  Under  these  assumptions  modal 
response  of  the  embedded  cylinder  can  be  described  by  one  generalized  coordi- 
nate, qn(t),  as  (Reference  9), 

w = 9n(t)  cos  n0,  (4-la) 

qn(t) 

v = sin  n0,  (4-lb) 

n 

where  w is  the  radial  displacement,  positive  inward  and  v is  the  tangential 
displacement  positive  in  the  direction  of  increasing  0 (counterclockwise). 

At  t = 0 the  n*"*1  flexural  mode  is  impulsively  excited  by  a suitable 
force  distribution  varying  like  a delta  function  in  time.  The  response  is 
composed  of  damped  harmonic  oscillations  (damping  due  to  radiated  elastic 
energy)  and  rapidly  decaying  starting  transients  (not  examined  further).  We 
wish  to  find  the  damped  frequency  and  damping  ratio  of  the  modes  of  such  a 
system  as  a function  of  cylinder  and  continuum  properties. 

The  solution  is  found  by  writing  the  single-degree-of-freedom  equa- 
tion of  motion  for  the  cylinder  subject  to  reactive  forces  from  the  surround- 
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ing  continuum.  Ignoring  starting  transients  we  look  for  solutions  of  the 
form 

. . —loot 

q„(t)  = ue 
n n 

Substituting  into  the  governing  differential  equation  gives  the  frequency  equa- 
tion, a numerical  solution  of  which  yields  the  admissable  complex  frequencies, 
w.  The  imaginary  part  of  u)  is  due  to  components  of  the  reactive  force  on  the 
shell  in  phase  with  the  shell  velocity. 

The  governing  shell  equation  is  derived  via  the  Euler-Lagrange  formu- 
lation. The  kinetic  and  potential  anergies,  T and  V,  for  lower  modes  of  the 
vibrating  cylinder  are: 


1 ■ 2 'O'  ah  (’  + S’)  V ■ i V' 

v- 112(1)’  l>!  - K ■>»’• 


(4-2a) 


(4-2b) 


where  m^  and  k^  are  the  modal  mass  and  stiffness,  respectively  and  prime  (') 
quantities  refer  to  the  properties  of  the  shell.  Substituting  into  the  Euler- 
Lagrange  differential  equation, 


* (~\  + 3V  _ p 

dt  ^3qJ  3qn  Fn 


(4-3) 


and  dividing  bv  m yields 
n 


q + U)  q = — , 
n n n m 

n 


(4-4) 


where  F is  the  generalized  force  due  to  radiation,  and  U)  2 = k /m  is  the 
n n n n 

in  vacuo  natural  frequency  of  the  cylinder. 

To  evaluate  the  generalized  force  of  the  continuum  on  the  cylinder, 
the  stresses  are  calculated  on  the  surface  of  a cavity  subject  to  the  cylinder 
displacements  (see  Figure  4-2).  Outward  propagating  waves  from  a driven  cavity 
are  found  in  terms  of  displacement  potentials  $ and  'f  (Reference  10)  as 


e 1 ‘ r A li  (K  r ) cos  n9  ; 
n n p 


e F)t  B H (K  r)  sin  n6  ; 
n n s 


K 


w/c 


K ; (o/c 
s s 


(4-5a) 

(4-5b) 


6 


where  A and  B are  complex  constants,  H is  the  Hankel  function  of  the  first 
n n n 

kind,  K , K are  wave  numbers  related  to  the  frequency  as  shown,  and  c and  c 
p s p s 

are  the  dilatational  and  shear  wave  speeds,  respectively.  Substituting  these 
solutions  into  the  displacement-potential  relations. 


u - $ + - S'  , 

r r r 0 


(4-6a) 


- - *P  , 

r 0 r 


(4-6b) 
-icot 


applying  the  prescribed  cavity  displacement,  Equation  (4-1)  with  q = u e 

n n 

yields  expression  for  and  as  functions  of  frequency.  Finally,  substi- 
tuting the  potentials  into  the  stress-potential  relations 


a = -A,  a + 2G'  (0  - -2  f 0), 

rr  c tt  rr  r 0 r ro 

P 


(4-7a) 


a 0 = 2G'  (~2  f - ¥ + - 

r0  2c  tt  rr  r 0r  r 

s 


* ) 


(4- 7b) 


results  in  the  normal  and  tangential  cavity  stresses  written  as 

a = 2G'  C (ft)  — e_it0t  cos  n9, 
rr  n a 


(4-8a) 


o Q = 2G ' D (ft)  — e ^t0t  sin  n0 , 
ru  n a 


(4- 8b) 


where  C and  D are  function  of  the  r.on-dimensionalized  frequency  defined  as 
n n 


ft  = &. 
c 

s 


(4-9) 


The  generalized  force,  F , is  found  by  the  method  of  virtual  work.  Referring 

n 

to  Figure  4-2  we  see  that 


2tt 

6w  = F <5q  = / (-(7  6 + 0 a 6 )ad0 , 

n n q rr  w r9  v 


(4-10) 


hence 


F = 2ttG'  (-C  (ft)  + — D (ft)  ) q 


n n 


(4-11) 


Substituting  F^  into  the  shell  equation.  Equation  (4-4),  and  setting 

q (t)  = u e iWt  yields  the  frequency  equation  for  the  cylinder  as 
n n 
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(4-12) 


3 2 + ~ ~7-nx  (-C  (ft)  + - D (12))  = 0, 

n Mn+1  n nn 


where  ft  is  the  damped  natural  frequency,  fi^  is  the  in  vacuo  natural  frequency, 
given  by 


_ 2 , . 2 

a co 


12a2 


n2  (n2  - 1): 
n^  + 1 


(4-13) 


with 


p'(l-G^)* 


and 


M.  2hP' 

ap 


(4-14) 


(4-15) 


is  the  mass  ratio  defined  as  the  ratio  of  cylinder  to  displaced  soil  mass  per 
unit  length. 

Writing  solutions  of  the  frequency  equation  as 

ft  = ft,  - iB  or  co  = co,  - i$, 
d d 

then  the  generalized  coordinate,  q^,  becomes 


q (t)  = u e 
n n 


-Bt 


(cos  co.t  + i sin  oj,t). 
d d 


(4-16) 


Damping  is  best  measured  by  the  damping  ratio,  q,  defined  as  the  ratio  of 

^d 

actual  to  critical  damping  (when  — r-  -*■  0)  and  given  by 

p 


D = sin  (tan  ^ — ) = - sin  (ARG(co)) 


(4-17) 


The  change  in  frequency  due  to  embedment  is  measured  by  t,,  the  ratio  of  damped 
frequency  to  in  vacuo  natural  frequency 

C = tOj/u  . (4-18) 

d n 

The  dynamic  response  is  explicitly  dependent  on  two-dimensionless 
parameters,  namely  the  in  vacuo  natural  frequency,  and  the  mass  ratio,  M. 

The  rigid  body  mode,  n = 1,  depends  on  the  mass  ratio  only,  while  the  flexural 
modes,  n = 2,  3,  . . .,  depend  on  both.  Observe  that  in  the  limiting  case 
of  a light,  flexible  shell  multiplying  the  frequency  equations  by  M and  letting 
M -*•  0 results  in 

-C  (ft)  + - D (ft)  = 0. 
n n n 
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The  solution,  is  the  frequency  of  an  oscillating,  inextensible  cavity. 

4.2.  SOLUTION  AND  APPLICATION  TO  CASES  OF  PRACTICAL  INTEREST 

Equation  (4-12)  was  solved  for  the  embedded  frequency  as  a function 
of  the  in  vacuo  natural  frequency  (Equation  4-13)  and  the  mass  ratio  (Equa- 
tion 4-15).  From  the  embedded  frequency,  the  frequency  (£)  and  damping  ratios 
(q)  were  calculated  (Equations  4-17  and  4-18).  Mode  1 (rigid  body),  Mode  2 
(ovaling)  and  Mode  3 were  investigated.  Except  for  very  massive  structures  the 
n = 1 mode  is  critically  damped,  while  the  n = 3 mode  is  qualitatively  similar 
to  n = 2 but  with  generally  less  damping;  these  modes  are  not  discussed  further. 
Typical  results  are  shown  graphically  for  the  n = 2 mode  with  mass  ratios 

M = 0.1,  0.34  and  1.0.  The  ratio  of  c to  c is  2 for  these  results.  The  n = 2 

P s 

mode  is  chosen  because  it  is  the  principal  mode  of  deformation  at  early  times 
for  many  protective  structures.  The  mass  ratios  cover  the  range  for  practical 
design  of  protective  structures.  For  example,  the  mass  ratios  for  the  arch, 
Reference  1,  and  for  rectangular  structure  3B  (£/h  = 10)  are  about  M = 0.4; 
rectangular  structure  3D  (the  subject  of  Section  2 of  the  present  report)  has 
mass  ratio  M - 1.0. 

The  solution  to  Equation  4-13  has  multiple  roots  for  each  mode  shape. 

This  arises  from  the  response  of  the  cavity  (not  of  the  shell).  The  consequence 
of  these  multiple  roots  is  multiple  frequencies  and  damping  ratios  of  the  coupled 
soil-structure  system.  The  root  which  dominates  the  response  in  each  mode  in  a 
particular  application  depends  on  the  nature  of  the  excitation  and  the  duration 
of  response  which  is  of  interest.  The  roots  selected  to  be  shown  in  this  report 
are  those  corresponding  to  results  reported  in  Reference  1 and  to  the  results 
obtained  by  the  present  writers  in  analyzing  the  box.  The  roots  may  not  always 
dominate. 

Figure  4-3a,b  show  the  frequency  and  damping  ratios  for  the  case  of  the 
ovaling  mode  with  mass  ratios  M = 0.34  and  1.0,  respectively.  The  abscissa  is 
the  frequency  parameter, 

au) 

Q = — . 
n c 

s 

Increasing  the  frequency  parameter  implies  larger  structure  (larger  radius,  a), 
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stiffer  structure  (larger  in  vacuo  frequency,  w ) or  softer  soil  (smaller  shear 

n 

wave  speed,  cg).  As  increases  the  damping  decreases  while  the  frequency  ap- 
proaches the  in  vacuo  value.  One  interpretation  of  this  is  that,  as  the  soil  be- 
comes softer,  the  soil-structure  system  tends  toward  the  in  vacuo  case. 

Comparison  is  made  in  Figure  4-3a  with  results  reported  in  Reference  1 
which  were  obtained  by  the  Doubly  Asymptotic  Approximation  for  the  same  soil- 
structure  system.  Agreement  is  satisfactory.  Comparison  is  made  in  Figure  4-3b 
with  results  of  the  present  finite  element  analysis  of  rectangular  structure  3D. 
This  comparison  is  only  approximate  because  of  differences  in  the  shapes  of  the 
structures  considered  in  the  two  analyses;  however,  ovaling-type  mode  shapes  are 
exhibited  by  both  structures.  Also,  in  the  case  of  structure  3D  there  is  a free 
surface,  which  tends  to  lower  the  embedded  frequency  and  decrease  radiation  damp- 
ing. In  view  of  these  differences  the  degree  of  agreement  is  acceptable. 

Figure  4-4  shows  the  damping  ratios  for  the  range  of  mass  factors  consi- 
dered; Figure  4-5  shows  the  corresponding  frequency  ratios.  This  is  the  same 
information  which  is  presented  in  Figure  4-3  with  the  case  M = 0.1  added.  These 
figures  can  be  used  in  the  design  and  interpretation  of  tests  and  calculations 
for  structures  whose  deformations  can  be  approximated  by  the  ovaling  mode  of 
a cylinder.  Two  practical  and  significant  applications  are: 

a.  A structure  and  its  backfill  can  be  designed  to  givg  maximum  or 

minimum  response  by  appropriate  choice  of  the  parameter  = — . 

s 

b.  Elastic  frequency  and  damping  ratios  can  be  estimated  for  single- 
degree-of-freedom  models. 
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SECTION  5 

SUMMARY  AND  RECOMMENDATIONS 

Modal  analysis  of  covered  and  uncovered  finite  element  models  of 
structure  3D  (£/ h = 4)  was  performed.  Due  to  accession  of  mass,  the  calculated 
natural  frequencies  of  the  covered  structure  were  less  than  frequencies  of  cor- 
responding modes  of  the  uncovered  structure.  Tests  on  the  physical  structure 
indicated  that  the  frequencies  were  essentially  unchanged  when  backfill  was 
added.  The  difference  between  measurements  and  the  computations  is  apparently 
due  to  arching  in  the  sand  backfill,  which  is  not  properly  accounted  for  in  the 
elastic  modal  analysis. 

Transient  analysis  of  a covered  model  of  structure  3D  was  performed  in 
which  two  different  backfill  wave  speeds  were  assumed  in  order  to  investigate 
their  effect  on  radiation  damping.  It  was  found  that  for  the  first  mode  the 
average  radiation  damping  for  the  entire  structure  can  be  expressed  as  15  per- 
cent of  critical  for  the  softest  backfill  considered  (cg  = 311  fps)  and  22  per- 
cent for  the  stiffest  backfill  (c  = 622  fps).  This  is  about  the  same  as  was 

s r 

found  in  a previous  analysis  of  an  arch,  in  which  damping  ranged  from  12  to  20 
percent  of  critical.  Radiation  damping  contributes  significantly  to  suppressing 
resonance  in  forced  vibration  tests  on  the  arch;  however,  resonance  in  the  rec- 
tangular structures  cannot  be  attributed  to  light  radiation  damping.  It  appears 
that  another  mechanism  such  as  arching,  possibly  peculiar  to  the  soil  type  used 
as  backfill  for  the  rectangular  structures,  plays  an  important  role  in  allowing 
resonance  to  develop  under  forced  vibration. 

Nonlinear  SDOF  models  of  the  rectangular  structures  were  developed  and 
subjected  to  loading  measured  during  a HEST  test  on  structure  3B.  The  loading 
included  the  HEST  cavity  pressure  and  the  soil-structure  interface  pressure. 

The  range  of  radiation  damping  ratios  obtained  during  the  transient  vibration 
studies  were  used  in  the  SDOF  models  to  investigate  the  effect  of  including  damp- 
ing on  the  response.  The  response  of  the  models  with  10  percent  critical  damp- 
ing were  30  percent  to  50  percent  less  than  those  of  the  undamped  models.  High- 
er damping  appears  to  be  justified  by  the  present  calculations  and  would  produce 
even  smaller  responses.  These  findings  strongly  indicate  the  need  to  include 
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radiation  damping  in  SDOF  models  of  shallow-buried  structures  in  which  deflection 
of  the  roof  results  in  deflections  of  the  walls  and  floor.  This  seems  to  be  the 
principal  mechanism  by  which  energy  of  vibration  is  transmitted  to  the  soil  and 
occurs  during  monotonically  increasing  deflections  of  the  structure  as  well  as 
for  vibrations. 

In  order  to  gain  insight  into  the  natural  frequencies  and  damping  charac- 
teristics of  embedded  protective  structures,  the  interaction  equations  for  the 
case  of  an  embedded  cylindrical  shell  were  solved.  The  results  are  presented  in 
graphical  form  such  that  the  ratio  of  embedded  in  vacuo  frequencies  and  fractions 
of  critical  damping  can  be  estimated  from  knowledge  of  the  in  vacuo  properties 
of  the  cylinder  and  the  properties  of  the  soil.  It  is  shown  that  these  results 
can  be  used  to  interpret  the  response  of  structural  shapes  other  than  cylinders 
(in  this  case,  a rectangular  structure)  providing  that  the  principal  mode  of  de- 
formation is  similar  (in  this  case,  the  ovaling  mode). 

One  remaining  uncertainty  is  whether  it  is  correct  to  assume  a constant 
fraction  of  critical  viscous  damping  based  on  stiffness  in  the  elastic-plastic 
regime  or  whether,  instead,  the  damping  should  vary  with  displacement.  Further 
research  should  be  aimed  at  this  problem  and  at  developing  simple  techniques 
for  estimating  the  radiation  damping  in  a wide  range  of  structural  configurations, 
depths  of  burial  and  backfill  wave  speeds  for  structures  of  interest. 
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APPENDIX  I 


STUDY  OF  DISCRETIZATION  USED  FOR  STRUCTURE 

The  finite  element  model  of  the  rectangular  structure  was  generated 
using  three-dimensional  continuum  elements  because  the  walls  and  roof  are  too 
thick  to  be  represented  by  thin  plate  elements.  A brief  study  was  performed 
to  determine  whether  the  lower  deformational  modes  of  the  roof  slab  could  ade- 
quately be  represented  by  a finite  element  model  having  two  continuum  elements 
through  the  thickness  or  if  three  or  more  elements  are  needed.  The  geometry 
of  the  test  problem  is  shown  in  Figure  1-1.  It  consists  of  a 6 feet  by  16  feet 
slab,  11  inches  thick,  which  is  simply  supported  on  all  sides.  By  applying 
appropriate  boundary  conditions,  only  one-quarter  of  the  slab  needs  to  be  con- 
sidered. 

Normal  modes  for  finite  element  models  of  this  slab  with  two  and  three 
elements  through  the  thickness  were  obtained  and  compared  with  each  other  and 
with  anlytic  expressions  for  thin  plates.  Typical  results  from  the  finite  ele- 
ment models  are  shown  in  Figures  1-2  and  1-3  where  the  frequencies  of  the  first 
and  second  modes  are  almost  exactly  equal  to  each  other.  Agreement  with  mea- 
sured frequencies.  Table  1-1,  is  difficult  to  evaluate  because  the  roof  of  the 
rectangular  structure  is  not  simply  supported.  Agreement  with  the  analytic  ex- 
pression for  frequencies  of  thin  plates 


tt2  / M2  , N2  \ 

[h.2  hj  V* 


where 


= 16  feet 


&s  = 6 feet 


k = flexural  stiffness,  ^ (1-v2) 


V = Poisson's  Ratio 


h = thickness  of  plate 


y = mass  per  unit  area 
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N = 1;  M = 1,  3,  5 for  modes  1,  2,  3,  respectively 

is  satisfactory  for  the  first  mode,  where  the  finite  element  model  has  a slight- 
ly lower  frequency.  This  is  because  the  finite  element  model  includes  the  ef- 
fects of  shear  deformation,  which  increases  the  flexibility  of  the  system.  For 
the  second  mode,  the  finite  element  model  has  the  higher  stiffness.  The  simplest 
explanation  of  this  finding  is  that  the  finite  element  model  is  too  stiff  to 
capture  the  second  mode  correctly  due  to  discretization  error  in  the  X-  and  Y- 
directions  as  well  as  through  the  thickness. 

The  results  of  this  study  indicate  that  the  analysis  of  the  rectangular 
is  just  as  valid  with  two  elements  through  the  thickness  of  the  roof,  walls  and 
floor  as  with  three  elements.  The  two  element  model  was  therefore  adopted  on 
grounds  of  economy. 


ss : simply  supported 

Thickness:  11  inches 


Figure  1-1.  Geometry  of  Slab  Used  to  Determine  Optimum  Discretization 
for  Roof . 
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TABLE  1-1 


r 


Mode 

1 

2 


FREQUENCIES  OF  ROOF  SLAB  (Hz) 

Finite  Element 
(Roof  Slab) 

147 

267 


Analytic 
(Roof  Slab) 

157 

196 


Experiment  (entire  structure)  237  Hz,  415  Hz,  470  Hz,  521  Hz 
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